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We have implemented three approaches to describe the thermodynamic properties of ferrimag- 
netic (S = 5/2, s — 2) spin chains. The application of cumulant expansion has been generalized to 
the ferrimagnetic chain in the presence of an external magnetic field. Using cumulants, we have ob- 
tained the field dependent effective Hamiltonian in terms of the classical variables up to the second 
order of quantum corrections. Thermodynamic functions, the internal energy, the specific heat and 
the magnetic susceptibility are obtained from the effective Hamiltonian. We have also examined 
the modified spin wave theory to derive the same physical properties. Finally, we have studied 
our model using quantum Monte Carlo simulation to obtain accurate results. The comparison of 
the above results and also the high temperature series expansion shows that cumulant expansion 
gives good results for moderate and high temperature regions while the modified spin wave theory 
is good for low temperatures. Moreover, the convergence regions of the cumulant expansion and 
the modified spin wave theory overlap each other which propose these two as a set of complement 
methods to get the thermodynamic properties of spin models. 

PACS numbers: 75.40.-s, 75.10.Hk, 75.10.Jm 



I. INTRODUCTION 



It will become evident that there are numerous highly 
interesting experimental systems which are effectively 
one-dimensional (ID) models. The ID models are more 
interesting from the theoretical point of view. The quan- 
tum effects which are highlighted in the ID spin models 
represent novel physical behavior. In this class, quan- 
tum ferrimagnets are mixed spin systems with antifcr- 
romagnetic interactions. Mostly, they are composed of 
the two type of spins, S ^ s. Two families of ferrimag- 
netic chains are described by ACu(pba)(H20)3 ■ T1II2O 
and ACu(pbaOH)(H 2 0) 3 .nH 2 0, where pba = 1,3- 
propylenebis(Oxamato) , pbaOH = 2 — hydroxo — 1,3 — 
propylenebis(Oxamato) and A = Ni, Fe, Co and Mn*£. 
Ferrimagnets, which occur rather frequently in nature, 
are somehow between the antiferromagnets and the fer- 
romagnets. Despite the fact that the homogeneous inte- 
ger spin chains show the Haldane gap in their low energy 
spectrum and the half-integer ones are gapless^, ID fer- 
rimagnets behave differently. The lowest energy band of 
the ID ferrimagnets is gapless which shows a ferromag- 
netic behavior while there is a finite gap to the next band 
above it which has the antiferromagnetic properties£*£i 
It is the acoustical and optical nature of excitations 
which is the result of two different type of spins in each 
unit cell. This behavior has been observed in the low 
and high temperature regime of quantum ferrimagnets 6 . 
There are many approaches to study the properties of 
the ferrimagnetic chains; The dual features of ferrimag- 
netic excitations can be illuminated by using the density- 
matrix renormalization group (DMRG)^ and quantum 
Monte Carlo (QMC) methods^. Numerical diagonaliza- 
tion, combined with Lanczos algorithmic and the scaling 



technique^, further have been applied to study the mod- 
ern topics such as the phase transition 1012 and the quan- 
tized magnetization plateau. 

The discovery of both ferromagnetic gapless and anti- 
ferromagnetic gapped excitations have led to the investi- 
gation of the thermodynamic properties. It has been pre- 
dicted that the specific heat at high temperatures should 
behave like an antiferromagnel£ which exhibit a Schot- 
tky peak at the intermediate temperatures. The modi- 
fied spin wave theory (MSWT) and QMC can be used to 
see this behavior however, QMC is not able to reach low 
temperatures sufficiently, to completely demonstrate the 
ferromagnetic behavior—. 

Most of the mentioned techniques such as QMC and 
DMRG have been used for ferrimagnets with small spins. 
The Hilbert space of large spins are growing exponen- 
tially and makes the computations more difficult. In 
recent years there are considerable attempts which has 
been focusing on the properties of new magnetic materi- 
als, such as magnetic molecules with large effective spins 
(S, s), or intermetallic compounds containing magnetic 
layers or chains. Using MSWT 13 and high temperature 
series expansion (HTSE)i4, one can describe the low tem- 
perature and high temperature properties of this systems, 
respectively. Moreover, the HTSE is not enough accurate 
even by taking into account higher terms (11th and 7th 
term for specific heat and susceptibility, respectively) . In 
addition the validity regime of HTSE is too far from the 
low temperature regime of MSWT to cover the full range 
of temperature. The mid temperature behavior of ferri- 
magnets with large values of S and s, has not received 
sufficient attention. In this respect, we used the MSWT 
to describe the thermodynamic behavior of the ferrimag- 
netic large spin chains. We have also implemented the 
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QMC simulation 15 as an accurate result for compari- 
son . At moderate temperatures, i.e Js < T < JSs 
(J is the exchange coupling) the results of MSWT do 
not coincide with QMC ones. Therefore, to describe 
the physical properties at mid temperatures we have em- 
ployed the cumulant expansion (CE) 1 ^ 1 T , Recently, this 
method has been used to study the finite temperature be- 
havior of large spin ferromagnetic and antiferromagnetic 
systema 18 ! 19 ! 20 ! 21 . 

In this article we have generalized the application of 
CE to obtain the magnetization and magnetic suscepti- 
bility of ferrimagnets. Moreover, the QMC simulation 
for (S = 5/2, s = 2) ferrimagnetic chain has been done 
to see the accurate behavior. We have observed good 
agreement between the QMC and CE results in the in- 
termediate and high temperature region. The outline of 
this paper is as follows: In Sec[H]we have employed three 
theoretical approaches, such as CE, QMC simulation and 
MSWT. The effective Hamiltonian and magnetization in 
the presence of the magnetic field are obtained using the 
CE up to the second order of (— ). In Sec lIIII the results 
and discussions have been demonstrated. We have com- 
pared our approaches with the high temperature series 
expansion (HTSE). It has been observed that CE and 
MSWT are two complement methods to get a good de- 
scription of large spin ferrimagnetic chains for the whole 
range of temperatures. 



II. THEORETICAL APPROACHES 

A. Cumulant Expansion 

Thermodynamic functions of any quantum spin system 
with a Hamiltonian H can be obtained by differentiation 
of the quantum partition function (Z) or its logarithm 
with respect to the appropriate parameters. Using the 
basis of spin coherent states Hn.i}) 2 ^, the trace of an 
operator is reduced to the integral over a set of classi- 
cal vectors. So the partition function is reduced to that 
of an effective classical spin system with the Hamilton 
function ^ 18 i 19 i 20 i 21 . The effective Hamiltonian TL can 
be expanded in terms of cumulantsi^ii! of the powers of 
H as follows, 

= 13(H {0) +H W +H i2) + ■■•)■ (2.1) 

where (3 — ^ and (0) c represents the cumulant of op- 
erator Gi§ii£. The function Tt evidently depends on the 
temperature, thus the calculation of the physical quanti- 
ties should be done with care. 

Let us consider the Hamiltonian (H ) of a ferrimagnetic 
chain which is composed of two kinds of spins S and s 



(S > s), alternatively. 
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where Hi's arc the external magnetic field at each sites. 
Expressing the spin operator on each site in the coor- 
dinate system with the z axis along the coherent state 
vector nf — rii 

n = n m +n (1) + n {2) 

+n^+n^+n {2 \ (2.3) 

where is the pure classical contribution, TiP^ and 
liP^ are the quantum corrections in the absence of a 
magnetic field (see Refs. |2£JEI). The field-dependent 
terms will be expressed in the following forms, 

JV JV 

i i 
jv jv 

{3lo 2 ' 2 

n h ] = 2«E^-w[( hai - 1 ' 11 «) 

id 

- (h 2 j_i • n 2l _i)(n 2i _i • n 2j ) + (h 2j • n 2l _i) 

- (h 2j • n 2i )(n 2 i-i • na,0] " ^ M 1 ~ 

i 

(h 2l _! ■ n 2l _!) 2 ) + 1 - (h 2l • n 2l ) 2 ] , (2.4) 

where u> = — . In the above expressions J = Js 2 and 
h = Hs are the exchange interaction and the reduced 
magnetic field, respectively. Again, TiP' is the classical 
contribution and the remaining higher orders are respon- 
sible for quantum corrections. The field-dependent terms 
of order h 3 , Jh 2 and h,J 2 in Ti.^ can be calculated with 
the help of cumulants corresponding to the mixed field- 
exchange terms. These terms are too lengthy and have 
been presented in the appendix. Quasiclassical expansion 
up to second order of (0(l/s 2 )) for the internal energy 
and the specific heat were investigated for different S, 
s in R.efs. |2f]ll2ll |. We will now calculate the CE of the 
magnetization and the susceptibility. In order to get the 
physical concept we have considered the nearest neigh- 
bor interaction and that the applied fields on each site 
are in the same direction, i.e hj = /in. Thus, the field- 
dependent terms of the effective Hamiltonian is reduced 
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to the following forms: 
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(2.5) 



where "J^ is expressed in terms of the classical vectors n 
and rii (see appendix). The partition function is repre- 
sented by the effective Hamiltonian defined in the previ- 
ous equations, 

(i - m { H ] + y^] 2 - (m {1) + f3 2 n^n^ 

+ y[W«] 2 - f3{H {2) + n { * ] ) + 0(l/s 3 ))] . (2.6) 

The reduced magnetic susceptibility is obtained by two 
times differentiating from the logarithm of the partition 
function with respect to h, i.e, 



dm 

X = hm — — . 

h^q dh 



(2.7) 



where H is the magnetic field and m is the scaled mag- 
netization per site and given by 



1 dlnZ 
Nd(Ph)' 



(2.8) 



In the limit of H — > 0, all of the terms containing h 3 or 
higher orders of h will vanish in the partition function. 
So, we will keep the expansion up to the second order of 
h in the partition function Ea. (|2.6|l . By the integration 
on the coherent states, we find the scaled magnetization 
in terms of the coupling and the magnetic field. The 
reduced magnetic susceptibility is as follow: 



X = 



P(1+lu 2 ) fl + B 
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(2.9) 



/3J and B = coth(w£) 



-K is the Langevin 



where £ 
function 

We have plotted in Fig. © the actual magnetic suscep- 
tibility (xa) of the (S = 5/2, s = 2) ferrimagnetic chain 
versus temperature. It is related to the susceptibility de- 
fined in Ea. H2.7| l by the following relation, 



,. 9m a 2 



(2.10) 



where m a = sm is the actual magnetization. The pre- 
sented CE result contains the quantum corrections up to 
the second order of (— ). We will discuss on the quality 
of our results in comparison with the other ones in the 
next section. 



B. Quantum Monte Carlo simulation 

We have implemented the quantum Monte Carlo 
(QMC) simulation for the ferrimagnetic (S = 5/2, s = 2) 
chain of length N — 64. The fairly large attainable 
length size gives us the thermodynamic properties as a 
very good approximate of the infinite size chain. In do- 
ing so, we have considered the Hamiltonian of Ea. (|2.2|) 
for N = 64 and without the magnetic field. We uti- 
lized the QMC algorithm based on the Suzuki- Trotter 
decomposition 23 of the Checkerboard- Type 2 ^. In this re- 
spect, we begin by breaking the Hamiltonian into four 

pieces, H = + H a + + Hi,, where H contains the 
interactions in z direction. H a and Hi, represent the in- 
teraction in the transverse direction alternatively. The 
partition function is expressed by the Suzuki- Trotter for- 
mula as the following, 



Z = Tre- pfl = lim Z„ 



fiHp fJH a f)H p6 b 

Z m : — Tr I e 2 >™ e "> e 2m e 



,(2.11) 



where m is Trotter number. Performing the trace oper- 
ation, we have a two-dimensional classical Hamiltonian 
rather than the one-dimensional quantum Hamiltonian. 
This classical Hamiltonian has 2m N spins. We have con- 
sidered the plaquette flip for the evolution of the Monte 
Carlo simulation. The reason is related to the huge num- 
ber of single spin flips which are not permitted because 
their Boltzman weight is zero. For instance, in the case 
of a plaquette of four spins which contains two S = | 
and two s = 2, there exists 900 different configurations. 
There are only 110 configurations with nonzero Boltz- 
man weights. All of these nonzero cases can be obtained 
by a plaquette ffipi^. The quantities like internal en- 
ergy, heat capacity, and magnetic susceptibility depend 
on the Trotter number (m). In the limit of m — ► oo, 
these quantities tend to their correct values. Therefore, 
It should be taken the biggest possible value of m, espe- 
cially at low temperatures. However, when temperature 
decreases the convergence relative to m becomes small. 
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In other words, using the big value for to makes two kinds 
of problems. Firstly, — becomes small so the state of 
the system (spin configuration) changes hardly (evolves 
slowly). Secondly, when m is large the global flips to 
change the total magnetization are hard to be accepted 
at low temperatures. Consequently, many Monte Carlo 
steps are needed to equilibrate the system at a large m 
and the low temperature. Therefore, for the mentioned 
reasons we have performed the calculations for each tem- 
perature with different values of m and utilizing the least- 
square extrapolation method fEq. ^2.12fl to find the limit 
of to — » oo) 

A(m) =A 00 + ^| + ^| + --- . (2.12) 

TO 171 

For moderate and low temperature regime (T < 3J) 
we have considered three different values for the Trot- 
ter number, to = 15, 20, 30. At higher temperatures, the 
convergence happens for the lower to values. To equlibri- 
ate the system we have spent 10 5 Monte Carlo steps and 
10 6 steps for measurement. Accuracy of the measured 
quantities depends on the temperature, for higher tem- 
peratures we got higher accuracy. However, the error bar 
is less than the symbol sizes in our plots. 

The internal energy, specific heat and magnetic sus- 
ceptibility of the ferrimagnetic (S = 5/2, s = 2) chain 
have been plotted in Figs. QJ, J5J) and © respectively. 
We will discuss our results in the next section when we 
compare it with the other ones. 

C. Modified Spin Wave Theory 

In the modified spin wave theory, usually it is consid- 
ered a single-component bosonic representation of each 
spin variable at the cost of the rotational symmetry. To 
simplify in the incoming calculations, we consider the 
following form of the Hamiltonian for the ferrimagnetic 
chain: 

N 

H=J^2(S i -s i - 1 +s i -S i ). (2.13) 

»=i 

Using the Holcstain-Primakoff and the Bogoliubov trans- 
formation, the Hamiltonian l|2.13|l is diagonalizcd Ti = 
-2NJSs + Ei + E a + Hi+H Q + 0(S'- 1 ), where E, gives 
the 0(S % ) quantum corrections to the ground state en- 
ergy and Hi is expressed in terms of a\ and and gives 
the quantum corrections to the dispersion relation, (see 
Ref.0])- a\ and f3' k are the creation operators of the 
ferromagnetic and antiferromagnetic spin waves with mo- 
mentum k, respectively. 

At finite temperatures, we assume that = 
n+ n ± Pj;(n",n + ), for the spin wave distribu- 
tion functions, where Pk(n~,n + ) is the probability 
of n~ ferromagnetic and n + antiferromagnetic spin 
waves appearing in the fc-momentum state and satisfies 



E„-,„+ Pk{n~,n+) = 1 for all fc's 13 . The substitutions 

= OiLo-k and = PkPk in the spin- wave Hamiltonian 
gives the zero field free energy, 

k 

+ E Pk{n- ,n+)\nP k {n- ,n+) . (2.14) 

k n — ,n+ 

To keep the number of bosons finite, one should apply 
the following constraint, 

(: S z ~s> :) = (S+s)N-(S+s)J2Yl ~ = °> ( 2 ' 15 ) 

j _ _i_ ^k 
k a—± 

where ujk = {(S — s) 2 + 4Ss sin 2 fc) 1//2 and the normal 
ordering is taken with respect to a and (3. By minimiza- 
tion of the free energy l|2.14|l with respect to Pk(n~ , n + )'s 
under the condition (|2.15|) we can obtain the free energy 
and the magnetic susceptibility at thermal equilibrium 
as follows 

F = E g + M (5 + s)N - J2 ln ( X + ( 2 ' 16 ) 

k <y=± 

k CT=± 

where n± = [ e (J"t-^s+s)/u h )/T _ ^ and { thc La _ 
grange multiplier to consider the constraint (|2.15l) . This 
set of equations has no closed analytic solution. In the 
case of (S = 5/2, s = 2) we have numerically solved Eqs. 
(|2.15|) and i|2.16|l in the thermodynamic limit, and visu- 
alized them in Figs. (JTJ, (J2J and ©. In previous equa- 
tions we have chosen ks — 1- and lu^ are the ferro- 
magnetic and antiferromagnetic excitation gaps, respec- 
tively. They have different values in the linear modified 
spin wave theory (LMSW) and perturbational interact- 
ing modified spin wave theory (PIMSW). In the PIMSW, 
the 0(5°) terms have been considered. Because the an- 
tiferromagnetic excitation gap is significantly improved 
by the inclusion of the 0(S°) correlation, the location of 
the Schottky peak can be also reproduced very well by 
the perturbational interacting modified spin waves. 

III. RESULTS AND DISSCUTIONS 

We have obtained the effective Hamiltonian of the ferri- 
magnetic chains in the presence of an external magnetic 
field to second order of cumulant expansion, Ea. l|2.5|l . 
The zeroth order term shows the classical contribution 
which simply represents the coupling energy of the clas- 
sical spins with the external magnetic field. Quantum 
corrections have a non-Heisenberg form and they are im- 
portant in the intermediate temperatures. 

In Fig.Q, we have shown the internal energy per unit 
cell of spins (nS-) versus temperature. The big difference 
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FIG. 1: Temperature dependence of the internal energy per 
cell of the ferrimagnetic (S = 5/2, s = 2) chain. The solid 
line: Cumulant expansion up to the second order, Dashed 
line: Classical part of the cumulant expansion, DashDot line: 
Perturbational interacting-modified-spin-wave-theory (PIM- 
SWT), Doted line: Linear-modified-spin- wave-theory (LM- 
SWT) and circle: Quantum Monte Carlo simulation with 
TV = 64 spins. 



FIG. 2: The specific heat per cell of the ferrimagnetic (S = 
5/2, s — 2) chain. Cumulant expansion up to second or- 
der (solid line), classical part of the cumulant expansion 
(dashed line), perturbational interacting-modified-spin- wave- 
theory (PIMSWT- dash-dot line), linear-modified-spin- wave- 
theory (LMSWT- doted line), high temperature series expan- 
sion (dash-dot-dot line) and quantum Monte Carlo simulation 
(open circle) with TV = 64 spins. 



between the zeroth order (classical contribution shown 
by dashed line) and the second order cumulant expan- 
sion (quantum corrections shown by solid line) shows the 
importance of the corrections in the intermediate and 
higher temperatures. The discrepancy is high, even in the 
present case of fairly large spins (S = 5/2, s = 2) which 
seems to behave classically. The reason is related to the 
dual features of ferrimagnets, i.c the low temperature be- 
havior is like ferromagnets and in the high temperatures 
they show antiferromagnetic behavior. There is a spec- 
tral gap in the subspace S to t — y (S — s) + 1, where the 
optical magnons play an important role. In the case of 
(S = 5/2, s = 2), the spectral gap of optical magnons at 
k = is; A = 1.36847J (Ref.|2gj). This means that the 
model does not behave pure classically. So, to describe 
the finite temperature behavior of the system, we should 
consider the quantum corrections to the classical part. 
At low temperature the second order of CE has a large 
deviation in comparison with the other results. Because 
in the low temperature region classical fluctuations are 
not strong enough to suppress the quantum ones. How- 
ever, the low temperature region has been excluded from 
the convergence domain by construction when CE is ex- 
pressed as a series in the order of (3Js < 1. Obviously, 
the classical term is dominant at very high temperatures. 

To have an impression on the accuracy of our results, 
we have plotted the results of the QMC simulation for 
comparison. The second order CE of the internal energy 
in Fig.(IJ fits very well on the QMC results for T > 2 J. 
This is actually the validity regime of our CE approach, 
T > Js. We have also plotted in Fig.Q the results of 
two different modified spin wave theory which deviate 



slightly from the QMC ones . However, the accuracy 
of the different schemes can be best visualized in the 
physical quantities like the specific heat and the magnetic 
susceptibility which is shown in Figs. ((21 El> ■ 

In Fig. J2J), we have plotted the specific heat per unit 
cell of spins for a (S — 5/2, s — 2) ferrimagnetic 

chain. The results of CE have been shown as the pure 
classical contribution and also the whole contribution to 
the second order. The big difference between them ver- 
ifies the significant corrections of the second order CE. 
We have also plotted the result of QMC simulation, for 
comparison. We observe very good agreement between 
the CE and the QMC results. According to the results 
presented in Figs. Q1I2JI, the QMC simulation results con- 
firm that the CE is a very good analytical approach to 
describe the thermodynamic properties of a ferrimagnetic 
system with large spins at moderate and high tempera- 
tures. 

We have also shown in Fig.J2J) the results of the MSWT 
for the specific heat. We have examined the LMSW 
and the PIMSW for our system. As it is observed 
from Fig. ©, both of LMSW and PIMSW can re- 
produce the high temperature behavior of heat capac- 
ity close to the QMC simulation results. Furthermore, 
they can show the Schottky peak at mid temperatures. 
Although the PIMSW can reproduce the location of 
the Schottky peak fairly well, it can not estimate the 
peak value well for large spins in comparison with the 
CE. The reason of this discrepancy in the MSWT is 
as follows. Let us come back to the spin wave theory 
and draw your attention to the bosonic Hamiltonian; 
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n 



-2NJSs + E 1 + E + Hi + Ho + CKS" 1 ) where 



k 



+J l (k){a k f3 k + al0l)], 



(3.1) 



and 7j(fc)'s have been introduced in Ref.[Tj|. The last 
two terms in Hi are the normal-ordered quassiparticlc 
interactions. In MSWT, whether LMSW or PIMSW we 
have eliminated these interactions, i.e we choose 



7l(*0 = O 



tanh 29 k = 



2v5scosfc 
S + s ' 



However, in the low and moderate temperatures, the 
magnon-magnon interactions play an important role. 
Therefore it is surmised that if we consider at least the 
first order of the quasiparticle interaction (i.e Jt(k) ^ 
and 7o(fc) = ), we can produce the Schottky peak value 
more precisely. Although the agreement between MSWT 
(without the quassiparticle interaction) and the other re- 
sults is not perfect, it is a remarkable success for the spin 
wave theory in the one-dimensional large spin ferrimag- 
net that all relevant features are quantitatively rather 
well reproduced over a very large temperature range. 

Finally, we have plotted in Fig. the magnetic suscep- 
tibility per total number of spins of the (S — 5/2, s — 
2) ferrimagnetic chain versus temperature. The result of 
the second order CE is shown by the solid line. The CE 
result shows qualitatively the features of the ferrimag- 
netic chains, i.e the quasiclassical x shows the divergence 
for T — > like a ferromagnet and a Curie law (~) de- 
cay at high temperatures. Meanwhile in Fig. J3J we have 
plotted the result of the QMC simulation. Our simula- 
tion result shows very well the antiferromagnetic feature 
of the ferrimagnetic system. But it can not produce the 
ferromagnetic behavior at low temperatures. The reason 
is as follow. The antiferromagnetic trait of the model 
does not depend on the size of system, i.e all antiferro- 
magnetic features can be reproduced in a rather short 
system size, while the ferromagnetic ones completely de- 
pend on the number of spins. This means that the fer- 
romagnetic features emerge only slowly with the growing 
of system size. In the QMC approach we have considered 
32 cells (64 spins) for simulation. The computation time 
grows exponentially by going to larger sizes. Specially, 
it happens in the calculation of magnetic susceptibility 
to reach the equilibrium condition. However, our main 
interest in this study is the results for intermediate and 
large temperature regions where reasonable values exist. 

We have also shown in Figs. J5J) and © the re- 
sults of HTSE for the specific heat and the suscepti- 
bility of (S = 5/2, s = 2) ferrimagnetic chain, respec- 
tively. The HTSE is an expansion in powers of (3J . 
Recently, Fukushima and his collaborators have imple- 
mented a suitable Pade approximation to obtain the ther- 
modynamic functions of the mixed spin chains. They 
have found the specific heat and the susceptibility up to 
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FIG. 3: Magnetic susceptibility per spin of the ferrimag- 
netic (S = 5/2, s = 2) chain versus temperature. The 
solid line: Cumulant expansion up to second order, Dashed 
line: Classical part of the cumulant expansion, DashDot line: 
Perturbational interacting-modified-spin-wave-theory (PIM- 
SWT), Doted line: Linear-modified-spin- wave-theory (LM- 
SWT), DashDotDot: High temperature series expansion and 
circle: Quantum Monte Carlo simulation with N = 64 spins. 



(0(/3J) n ) and (0(/3) 7 ), respectivljiii The HTSE results 
for the specific heat shown in Fig.J2Jl converges to the 
QMC results at T > 4J. However, the deviation from 
the QMC result is more pronounced for the magnetic 
susceptibility shown in Fig.©. 

It is worth to mention the two differences between the 
CE and HTSE results. Firstly, the convergence region 
of CE is larger than the HTSE one, ie, the CE is valid 
for T > Js wh ereas the validity of the HTSE is for 
T > y/S(S+l)J. Secondly the HTSE fails to produce 
the Schottky peak of the specific heat While the CE can 
generate it as well as the QMC simulation. 

Our results state that the combined methods of the 
cumulant expansion for T > Js and the modified spin 
wave theory for T < Js, give a good approximation for 
the whole finite temperature behavior of quantum fer- 
rimagnets. The reason is related to the overlap of the 
convergence regions of the mentioned method. 
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APPENDIX A 

(2) 

The expressions of the second term (H h ) of the field 
dependent effective Hamiltonian are obtained as follows, 

N 
2 

*i = ^{wn-n 2l -i(l-(n-n 2 ,_i) 2 ) 
+n • n 2 i(l - (n ■ n 2i ) 2 )}, 

N 

1=1 

- rii ■ n j;+ i - (n ■ rij)(n • n i+ i) + (n,; • n i+ i) x 
((n. ni ) 2 + (n.n 4+1 ) 2 )}, 

JV 

f 3 = ^(n< n i ){(n i _i • n;) 2 + (n* • n i+1 ) 2 } 

i=i 

jv 

2 

+ ^ {2(n • n 2j -i)(n 2 ;_i • n 2i )(n 2 ,_i • n 2i _ 2 ) 

i=i 

+ (n • n 2i )(n 2l • n 2i _i)(n 2i _i • n 2i _ 2 ) 
+ (n • n 2i _ 2 )(n 2 j_ 2 • n 2i _i)(n 2i _i • n 2l ) 

- (n • n 2i _ 2 + n • n 2l -i + n • n 2i ) x 

+ (n • n 2i _ 2 )(n 2i _ 2 • n 24 _i)(n 2l _i • n 2i ) 

- (n • n 2i _ 2 + n • n 2i -i + n • n 2i ) x 
(n 2 j_i • n 2 i + n 2 i_ 2 • n 2 j_i) 

- (n 2i • n 2i _ 2 )(n- n 2i + n- n 2 j_ 2 )} 



* 4 = ^(n-n,)((n 4 _ 1 -n l ) 2 + (n,-n l+1 ) 2 ) 

i=i 

jv 

2 

+ ^2 {2(n • n 2i )(n 2i ■ n 2i _i)(n 2i ■ n 2i+1 ) 

i=l 

+ (n • n 2i _i)(n 2i _i • n 2l )(n 2l ■ n 2l+1 ) 
+ (n • n 2l+ i)(n 2l+ i • n 2l )(n 2l _i • n 2i ) 

- (n • n 2l _i + n • n 2i+1 + n • n 2i ) x 
(n 2l _i • n 2i + n 2i • n 2i+1 ) 

- (n 2i _i • n 2i+ i)(n • n 2i _i + n • n 2i+ i)} 

2 

*5 = [(l-n 2i _i ■n 2i )(2n-n 2l + 2n-n 2t _ 1 

i=l 

+ (1 - 3n 2l _i • n 2l )(n • n 2i _i + n • n 2i )) 
+ (1 - n 2j _i ■ n 2l _ 2 )(2n • n 2l _ 2 + 2n • n 2i _i 
+ (1 - 3n 2j _i • n 2l _ 2 )(n • n 2l _i + n ■ n 2l _ 2 ))] 
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